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Abstract 

The time- dependent Ginzburg-Landau formalism is applied to analyze the 
vortex states and vortex dynamics in superconducting spherical nanoshells, 
subjected to mutually perpendicular strong dc and weak ac magnetic fields. 
We demonstrate that nonuniformity of the shell thickness can dramatically 
affect the ac magnetic response of a 3D array of superconducting nanoshells. 
Remarkably, this response is strongly influenced not only by the relevant 
geometric and material parameters and the ac-field frequency but also by 
the magnitude of the applied dc field: by changing this field the real part 
of the effective ac magnetic permeability can be tuned from positive values 
significantly larger than one down to negative values. 
Keywords: 

A. Superconductors; A. Nanoshells; D. Vortices; D. Magnetic permeability 
1. Introduction 

Over the last decade, advances in nanotechnology have inspired experi- 
mental and theoretical studies of vortex behaviour in superconductors with 
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ferro- 



superconducting spheres 



15| . including those with a 



magnetic dipole inside(l6l. Il7j. Nanoshells are hybrid nanostructures, which 
typically consist of a dielectric core (usually a silicon oxide nanograin) , coated 
with a thin layer of metal {lil , [19I ] . When the metal in its bulk form is a su- 
perconductor, the nanoshell below the critical temperature will also exhibit 
superconductivity in the thin shell around the insulating core. It has been 
shown that those nanoshells allow the coexistence of a Meissner state and a 
vortex state in equilibrium on one and the same superconducting layer 15]. 



In superconducting spherical nanoshells the surface curvature and the applied 
homogeneous magnetic field lead to a Magnus-Lorentz force, which pushes 
the vortices and antivortices 20J towards the opposite poles of the shell. The 



same effect can be expected in the absence of any external magnetic field for 
spherical superconductor-ferromagnet nanoshells with uniformly magnetized 
core. This can be considered as an effective pinning of vortices and antivor- 
tices at the poles, which stron gly affects both the equilibrium distributions 



of vortices and their dynamics 
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In this Communication the vortex dynamics in superconducting nanoshells 
is investigated for the case when in addition to a (strong) dc magnetic field, 
which provides vortex pinning, also a weak ac magnetic field is present. In 
particular, we analyze the ac magnetic response of nanoshells with nonuni- 
form thickness of the superconducting shell. We show that the effective ac 
magnetic permeability of nanoshells, arranged in a 3D array, is very sensitive 
to the applied dc magnetic field. For one and the same array of nanoshells, 



2 



the real part of the ac magnetic permeability at a given frequency can be 
tuned from relatively large ("superparamagnetic") positive values at rela- 
tively high dc fields to negative values at lower dc fields. 



2. Model 

A sketch of a spherical nanoshell with thickness d and radius R is shown 
in Fig. Ufa). A vortex-antivortex pair is induced in the nanoshell by a ho- 
mogeneous dc magnetic field Bd c parallel to the z-axis. An additional oscil- 
lating magnetic field B ac = e x B UJ cos (ut) with frequency u> and amplitude 
B u <C -Bdc is applied along the x-axis. In our calculations we use spherical 
coordinates r, 9, <fr, where 9 = corresponds to the positive direction of the 
z-axis in the Cartesian co-ordinate frame. We consider a spatially periodic 
3D rectangular lattice of identical nanoshells, which are coupled to each other 
only electromagnetically. The unit cell of the lattice contains one nanoshell 
in the unit-cell center and has sizes L x , L y , L z > 2R, so that the edges of 
the unit cell do not intersect with the nanoshell. The amplitude of the ac 
field is kept relatively small (B u < 10~ 2 ) so that the response of nanoshells 
to this field is linear. 

To describe the vortex dynamics in nanoshells, we exploit the time- 



dependent Ginzburg-Landau (TDGL) approach. Like in 21[, the relevant 
quantities are made dimensionless by expressing lengths in units of \/2£, 
time in units of 7t/i/[4/cb(T c — T)] w 11.6tgl, magnetic field in units of 
$o/(47r£ 2 ) = H c2 /2, current density in units of $o/(2\/27r/xoA 2 £), and scalar 
potential in units of 2ks{T c — T)/ (7re). Here, $o = 7r^/e is the magnetic flux 
quantum, /x is the vacuum permeability, A is the penetration depth, £ is the 
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Figure 1: (Color online) (a) Sketch of a nanosliell with uniform thickness d. The red and 
blue shadow schematically show a vortex and an antivortex, respectively, (b) Distribution 
on the order parameter in a nanoshell with R = 2, d = 0.35, k = 0.35 for L x = L y = L z = 
5.2, i?dc = 0.6 and B ac = 0. (c) Distribution of the radial component of the normalized 
magnetic field in a nanoshell with R = 2, d = 0.35, k = 0.35 for L x — L y — L z — 5.2, 
B dc — 0.6, B u — 0.003, uj = 0.003. Panels (d) to (f): same as in panels (a) to (c), 
respectively, but for the case of the nonuniform shell thickness described by Eq. (fT2"]) with 
da = 0.35, p — 3 and q = 5. 
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coherence length, tql is the Ginzburg-Landau time, and H c2 is the second 
critical field. 

We consider the case where the thickness of a superconducting shell is 
significantly smaller than the coherence length £, so that variations of the 
order-parameter magnitude across the shell as well as currents in the radial 
direction are negligible. In general, the thickness d is assumed to be nonuni- 
form. In this case, using the aforedescribed units the TDGL equation for the 



order parameter if) [22| . adapted to the case of spherical geometry, takes the 
form 



d\Rsm9d(f) 9 J \Rsin9 

1/10 cote a 

where <p is the scalar potential and A is the vector potential. 

The scalar-potential distribution is determined from the condition 

V-j = 0, (2) 

which reflects the continuity of currents in a superconducting shell. The 
total current density j in the shell is given by the sum of the normal and 
superconducting components: 

j=jn+js, (3) 



a 



j, = Im«>*VV0-A|./f, (5) 



where a is the normal-state conductivity, which is taken as a = 1/12 in our 



units (see Ref. 



23(|) . The vector potential A, for which we choose the gauge 



V • A = 0, can be represented as 

A = A dc + A ac + A s . (6) 

Here A^c denotes the vector potential corresponding to the dc magnetic field 
Bdc, while A ac is related to the ac magnetic field. The vector potential A s 
describes the magnetic fields, induced by the currents j, which flow in the 
superconducting shells: 

where k = A/£ is the Ginzburg-Landau parameter and integration is per- 
formed over the volume of the superconducting shells. The results, described 
below, are obtained for a fixed value of the Ginzburg-Landau parameter 
{k = 0.35). 

When solving numerically the TDGL equation (JTJ, the gauge invariance of 
the discretized expressions is preserved by introducing link variables following 



the method of Refs. 



24 



25| . Two-dimensional grids, used in our calculations, 
typically have > 30R equally spaced nodes in the ^-interval from to n and 
> 45.R equally spaced nodes in the ^-interval from to 2tt, so that the 
distance between the neighboring nodes does not exceed 0.2£ (recall that the 
dimensionless parameter R is given by the ratio of the nanoshell radius to 
\/2£)- The step h t of the time variable t is automatically adapted in the 
course of calculation. This adaptation is aimed to minimize the number of 
steps in t and - at the same time - to keep the solving procedure accurate: 
the value of h t is chosen as large as possible provided that for t new = t Q id + h t 



the results, calculated with the steps h t and h t /2, are practically the same. 
Typically, the step h t is ~ 1CT 5 to ~ 1CT 3 depending on a specific distribution 
of the order parameter. For momentary distributions of the order parameter 
and the vector potential A, an iteration procedure is used to determine from 
Eq. ([2]) the corresponding distribution of the scalar potential p with a relative 
accuracy not worse than 10 -4 . The (time-dependent) vector potential A 
and the corresponding link variables are calculated using Eq. (J7J), where the 
contributions of ~ 10 4 neighboring nanoshells are taken into account. 



According to Ref . 26] , the effective relative magnetic permeability, which 
corresponds to the response of a nanoshell lattice under consideration to the 
ac magnetic field B ac = e x B w cos cut, can be expressed as 

^ = R +/r( c )\ , ,7r( s )\ ' ^ 
+ \Dsw jx + 1\-Dsuj )x 

where the numerator originates from the ac -B x -field, averaged over a unit-cell 
face, normal to the direction of the applied ac field, while the denominator 
is determined by the total ac if^-field, averaged over a unit-cell edge parallel 
to the applied ac field, times fio- Here the quantities 

B S = ^ / dt'B s (t')cos(cut') (9) 
t 

and 

t+2n/uj 

Bi sJ = ^ J dt'B s (t')sm(cot'), (10) 



t 



with B s = V x A s , describe the in-phase and out-of-phase magnetic field, 
respectively, induced by the superconducting nanoshells at frequency u. As 



seen from Eq. (J7J), these fields are inversely proportional to k 2 . By introducing 
the normalized complex field with the real part = a^/B^ and the 
imaginary part b" = ~Bj£j/B w , equation (JSJ) simplifies into 

= a + (u.r 1 • (ii) 

3. Results and discussion 

In Fig. UXb), we plot the equilibrium distribution of the order parame- 
ter in a nanoshell with radius R = 2 and uniform thickness d = 0.35 at 
Bdc = 0.6 and _B ac = 0. This distribution corresponds to a single vortex- 
antivortex pair with the vortex (antivortex) pinned by the Magnus-Lorentz 
force to the northern (southern) pole of the nanoshell. In response to an ac 
magnetic field with B u = 0.003, u = 0.003, applied to a cubic lattice of those 
nanoshells with L x = L y = L z = 5.2, the nanoshells induce their own mag- 
netic field. The radial component of the corresponding normalized field b^ in 
a nanoshell is shown in Fig. HJc). In the equatorial regions of the nanoshell, 
the magnitude and the sign of the field (b'J) r are mainly determined by the 
diamagnetic Meissner currents, induced in the nanoshell by the applied field 
B ac . The distribution of (W u ) r near the poles is governed by the oscillatory 
motion of the vortex and antivortex. Indeed, the field B ac , parallel to the 
tends to shift the (anti)vortex from its equilibrium position on the 
z-axis towards the x-axis. Importantly, each (anti)vortex carries a relatively 
large magnetic moment, which is (anti)parallel to the r-axis. The local mag- 
netic fields created by a vortex (antivortex) in the nanoshell and its vicinity 
are comparable in magnitude to the relatively strong dc field B^ c . There- 
fore, oscillations of (anti)vortices can produce an appreciably high ac field. 
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Figure 2: (Color online) Distribution of the x-component of the normalized magnetic field 
along a unit-cell edge, parallel to the x-axis, for nanoshells of uniform thickness d — 0.35 
(solid line) and nonuniform thickness described by Eq. (fl2|) with do = 0.35, p = 3 and 
q = 5 (dashed line). The calculations correspond to R = 2, k = 0.35, L x = L y = L z = 5.2, 
B dc = 0.6, B u = 0.003, and u> = 0.003. 

However, due to strong pinning of (anti)vortices to the nanoshell poles, the 
oscillation amplitude of (anti)vortices in the nanoshell lattice under consider- 
ation is rather small. As a result, the normalized magnetic field on the 
unit-cell edge parallel to the x-axis (see the solid curve in Fig. E]) is actually 
dominated by the contribution of the Meissner currents. Consequently, the 
effective magnetic permeability of the nanoshell lattice /i^ ~ 0.915, deter- 
mined by Eq. (flTT) . corresponds to a diamagnetic medium. 

For nanoshells with a nonuniform thickness, the situation can become 
significantly different from that described above. Here we restrict ourselves 
to the case when the shell thickness d varies according to the formula 

f (l + cos p fl) forf5<7r/2 
d = d { (12) 

[ (1- |cos0|«) for 6 > tt/2 

with p, q > 2, as it is schematically shown in Fig. d^d). Those variations 
of d enhance the anivortex pinning at the southern pole of a nanoshell and, 
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at the same time, tend to expel the vortex from the northern pole. If these 
variations are sufficiently large, the equilibrium position of the vortex in an 
isolated nanoshell corresponds to a ring-like potential well at 9 ^ 0. Since this 
potential well has perfect symmetry, the vortex can easily move within 
the well (in the direction) when applying even a very weak additional 
force, caused by the field B ac . In a rectangular nanoshell lattice the 
symmetry of the potential well for a vortex is somewhat violated, mainly 
due to the magnetic field of the Meissner currents induced by the field B dc in 
the neighboring nanoshells. In such a lattice, subjected to a weak ac magnetic 
field parallel to the x-axis, the vortex typically oscillates around <fi = tt/2, 
as illustrated by Figs. [He) and (f), or around <ft — 3ti/2. At relatively low 
frequencies u the amplitude of the vortex oscillations can be rather large, and 
the radial ac magnetic field induced by these oscillations appears by orders 
of magnitude stronger than the applied field [see Fig. [TJf)]. The ac field, 
induced by the nanoshells on the unit-cell edge, parallel to the x-axis, is fully 
dominated by the contribution of the oscillating vortices. In the central part 
of the unit-cell edge the x-component of this field is opposite to the applied 
ac field and significantly exceeds it (|(^)x| > 1; see the dashed curve in 
Fig. |2j). so that the total normalized field 1 + (W u ) x becomes negative. 

In Fig. [3J a) the x-components of the normalized magnetic fields and 
b" averaged over the unit-cell edge, parallel to the plotted as a 

function of u. As seen from Fig. 13(a) , at u < 0.04 the averaged field {b' u ) x is 
negative. Moreover, for u < 0.004 its magnitude is larger than 1, so that the 
real part of the effective magnetic permeability of the nanoshell lattice /jf u , 
defined by Eq. ( fTTj) . appears negative [see Fig. [3(b)] . As further seen from 
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Fig. El in the low-frequency range the out-of-phase field component (tQ x and, 
consequently, the imaginary part of fi^ rapidly decrease when decreasing uj. 
This decrease originates from a reduction of the energy dissipation by the 
vortices, which move with smaller average velocities at lower frequencies. At 
relatively high frequencies, the amplitude of the vortex oscillations and their 
contribution to (bj) x and /i w diminish with increasing u. For u > 0.01 the 
magnetic response of the nanoshell lattice is dominated by the contribution 
of the dissipationless Meissner currents, induced by the field B ac . In this 
frequency range the effective permeability /i w corresponds to a diamagnetic 
medium, and the imaginary part of /j, w gradually vanishes with increasing uj. 

Interestingly, the effective ac magnetic permeability of a nanoshell array 
is sensitive not only to the relevant geometric and material parameters and 
the ac-field frequency but also to the magnitude of the applied dc field B& c . 
In this connection, it seems worth mentioning that in nanoshells with R ~ 2 
the state with one vortex- antivortex pair appears (meta) stable in a rather 
wide range of B& c [14]. As illustrated in Fig. HI the ac magnetic response 
of a nanoshell array at a fixed, relatively low frequency can be efficiently 
controlled by simply varying the dc field B& c within this range. Indeed, 
the plots in Fig. ID^a) show that with increasing B<\ c the magnitude of the 
negative normalized field (b'^jx changes from relatively large values | (p' U] ) x \ > 1 
to values < 1- This behavior can be explained by a shift of the 

equilibrium vortex position towards the nanoshell pole and the corresponding 
decrease of the angle between the vortex magnetic moment and the z-axis 
when applying a higher dc field B^ c . As seen from Fig. HJ^a), this leads to 
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Figure 3: (Color online) (a) x-components of the normalized magnetic fields and b" 
averaged over the unit-cell edge, parallel to the x-axis, as a function of the applied ac field 
frequency, (b) Real and imaginary parts of the effective magnetic permeability of the array 
of nanoshells as a function of the applied ac field frequency. The calculations correspond 
to R = 2, d = 0.35, p = 3, q = 5, k = 0.35, L x = L y = L z — 5.2, B^c = 0.8, and 
B^ = 0.001. 
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Figure 4: (Color online) Averaged x-components of the normalized magnetic fields and 
b" (a) and real and imaginary parts of the corresponding effective magnetic permeability 
fj, u (b) as a function of the applied dc magnetic field. The calculations are performed for 
L x = L y = L z = 5.2, p = 3, q = 5, to — 0.001 (case A) and L x = L y = 4.4, L z = 6, 
p = q = 5, u) = 0.003 (case B). Other relevant parameters are R — 2, do — 0.35, k = 0.35, 
B u = 0.001. 

an evolution of the ac magnetic permeability from negative values at 
lower dc fields B& c to relatively large positive values, which correspond to 
(super)paramagnetic ac response of the nanoshell lattice at higher fields B dc . 



4. Conclusions 

Using the time-dependent Ginzburg-Landau approach, we have investi- 
gated vortex states in superconducting spherical nanoshells with nonuniform 
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thickness, subjected to a homogeneous dc magnetic field, and analyzed the 
vortex dynamics in the presence of an additional weak ac magnetic field, per- 
pendicular to the dc field. It is shown that by increasing the shell thickness at 
the pole, determined by the direction of the dc field, the equilibrium vortex 
position can be shifted from this pole to the surrounding ring-like potential 
well. In this state, the amplitude of vortex oscillations in an ac magnetic field 
appears strongly enhanced, so that the magnetic response of a nanoshell is 
fully dominated by the contribution of the oscillating magnetic moment of 
the vortex. As distinct from the case of nanoshells with uniform thickness, 
which are typically characterized by a diamagnetic ac response, for a 3D 
lattice made of nanoshells with nonuniform thickness the real part of the ef- 
fective ac magnetic permeability can take negative values at sufficiently low 
frequencies. We have demonstrated that the ac magnetic response of such a 
lattice can be efficiently controlled by the magnitude of the applied dc mag- 
netic field. By varying this dc field the real part of the effective ac magnetic 
permeability of the nanoshell lattice can be tuned in a wide range: from 
negative values to relatively large positive ( "superparamagnetic" ) values. 
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